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Summary 


The paraxial trajectories in electron lenses are derived from an integral 
equation. The Liouville-Neumann expansions of the solutions of this 
equation Icad to expressions for the magnification, the focal distances and 
the positions of the focal points and cardinal points. The numbers of integra- 
tions to be performed in the individual terms of the expansions used to 
describe lens characteristics are reduced to half, as compared with the 
normal treatment. The focal and cardinal points are defined as osculating 
elements similar to those introduced by Glaser. 


§ 1. Introduction. The paraxial theory of rotationally symmetric 
clectron lenses is based on the neglect of second and higher powers 
of the distance 9 from the axis when deriving the cquation of motion 
of the’electrons. The linear differential equation thus obtained is of 
the following form for electrons originally in a meridional plane: 


(1) 


The coordinate z refers to distances in the direction of the symmc- 
try-axis. In the most gencral case of a combined electrostatic and 
magnetic lens the dependent variable is given by y = o(z) V# (0, 2), 
V (0, z) being the potential in e.s.u. of the electric field. The cocffi- 
cient #?(z) depends on the electrostatic potential V(z) and the 
magnictic field H(z) along the axis according to ') 


3 aay e FH? 
~ 16 V2 © Bme? V’ 


(2) 


(2) 


= #16 5 
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The effect of space charge is neglected here, 

Only in a few exceptional cases is k?(z) of such a form that (1) can 
be solved analytically for a ficld approximately realized in practice. 
In this connection we mention in particular the case 

ke 


*@) = 7s 


(1 + 27/4)? 
that corresponds to the bell shaped magnetic field 


which has been investigated thoroughly by W. Glaser?), This 
field leads to expressions for e(z) that depend on trigonometric 
functions, viz. : 


2\3 
ae (1 ae 5) sin {c + ha)! (arc tan — oh. , 


Further, a magnetic field corresponding to electron trajectories 
described by Legendre functions is given by 


Hy 
cosh (z/a) ° 


The most general solution in this case, viz, 


e= AP, {tanh ({)} + BQ, {tanh (=) : 


depends on Legendre functions, the degree of which follows from 
' the relation : 
ea’ FR? 
8me? V 

As a rule, however, it is necessary to resort to some numerical 
method or other in order to derive results from (1) when the field 
distribution along the axis is given by experimental data. In parti- 
cular, approximation methods introducing finite differences instead 
of differential quotients have been frequently applied. A special. 
method based on successive approximating integrations has been 
worked out by M. v. Ments and J.B. Le Poole’). In the 
present paper the méthod used by these authors is simplified by 


yy + 1) = 
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lowering the number of integrations occurring in the various correc- 
tion terms. The treatment is based on an integral equation discussed 
in the next sections. 


§ 2. Derivation of the paraxial integral equalion. I.ct us consider 
the general case of an immersion Iens. The valuc h?(z,) of 4?(z) at the 
object plane = = 2) may thus differ from zero. The equation to be 
derived will refer to a special clectron trajectory, characterized by 
the values vp and ygof y and dy/dz respectively at the object plane. 
According to the relation y = eV! these quantilics are proportional 
to the distance og from the axis and to the tangent of the angle with 
the z-axis at the object point from which the electron is starting. A 
first integration of the paraxial equation, viz. d?y/d2? = — h(z) ¥, 
then yields ‘ 

dy 
de 


2% Vo a | ds h(8) y(2). 


=a 
A second integration gives 


; 7 ray zs z 
v=o + ¥9 (2 — 20) —Sds, fds MS) ¥@). 
co 20 


By inverting the order of integration the double integral proves to 
be equal to 
fas PS) yO) fad, = Sas HS) yO) 4). 
ay 4 =0 
Il is to be noted that this reduction to a single integral halves the 
numbers of inlegrations occurring in the expansions lo be derived in the 
following sections. 
We thus arrive at the following integral equation: 


vl) = yo + yole— ae) —Sas BEV E—D yO, 


which completely determines the paraxial behaviour of the lens 
system, : 


§ 3. Solution of the integral equation by the Liowville-Neumann 
method. The first step in this method +) consists in substituting into 
(3) for y(é) its value according to the equation itself, that is 


v(C) = yo + 95 (6 — 20) — Sab, HUE) E— 8) ¥EE. 
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The evaluation of this first substitution yields 


yl2) = Yo + y3lz — 29) —S de #0) (@ —9 fy +.) (6 —20)} + 


F t 
+ fae h(E) (2 — 6) fad, (61) (6 — &) v(Q). (4) 


20 
The second step amounts to substituting for y(é,) in the last integral 
once again its value according to (3), that is 


& 
y(E1) = Vo + Yo (0) —- %) —f dé, RE) (2, — &2) y(C.). 


The third step consists of a further application of (3) to y(é,), and so 
on. By repeating this procedure ad infinitum we obtain the scries 


y(2) = ¥9 + 9h (2 — 2) —SaE PO) (@—9) fy +94 (C—%)} 4 


2 c 
+ f db RC) (2 — 0) fad, (6) (6 f1) {Yo + Yo (S1 — 20)} — 


z t ray 
— fat Wf) (2 — 0) fey WE) (6 — b1) f de, R(E,) (6, — &2) - 
Yo + Vo b2— Ao} +-.-- (8) 


This expansion with the aid of so-called ttevated kernels represents 
the exact solution of (3) in the case of convergence. A sufficient condi- 
tion for convergence is that the product of the integration interval 
and the maximum value of the main kernel should not exceed unity. 
We are only interested in the space bounded by the object plane 
2== 2, and the corresponding image plane z = z,;. The main kernel being 
given by 4?(¢) (e—2), with a maximum valuc max k?(z) x (2; —%)), 
the above sufficient condition for convergence amounts to 


max k*(z) x (2; — 2)? <1. 
Ifor a magnetic lens this reduces to . 


8 2 
re ae (2; Zo)? < ia V; 
é 


2 
i. max 


(2; — %)? < 48.2 V 


if Hise, %— % and V are expressed in gauss, centimetre and volt 
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respectively. The condition (6) can be given in a morc illustrative 
form by replacing the magnetic field 7 by the curvature 

,__ mue 

eH 
of the trajectory of an electron moving with velocity v perpendicular- 
ly to a homogencous field of the same strength. The cquality (6) 
then reads 

a = 2 < 2 Kt 


The expansions to be discussed in the next sections therefore con- 
verge more rapidly according as the quantity (z;— 2))/2 R,,,, iS 
smaller. 


§ 4. Particular solutions. Determinalion of the image plane. The 
expansion (5) disintegrates automatically into 


(2) = Vo Vi(2) + Yo Volz), (7) 


in which occur two particular solutions y,(z) and y,(z) of the original 
equation (1). These special solutions are given by 


2 2 4 
¥1(2) = 1 — Lae re) (z—6) ties #(S) Go) fade, REM E— 


— faz HE) e—9) fab, PE) CE) fae ED &—A) + --- 
and 


Yo{2) = 2 — 2 — 


— saree Ne De — a) + Jar Hl Neat, P(E) (C—C4) (6, -%0) — 


jaz h(E) Mes ag UC) )(S—é fatal (2) (€:—22) 62— 


The solution y,(z) = V#(z) @,(z) is characterized by 4,(z) = ! 
and ¥i(z) = 0. Therefore y,(z) represents the trajectory of the 
electrons leaving the object at a point, a distance g@,(%) = 1/ VE (2) 
away from the axis in a direction to be determined from y,(z) = 0; 
for magnetic lenses (V independent of z) this dircction is parallel to 
the axis. The function y,(z) = V}(z) 9,(z) likewise describes the 
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elbctrons leaving the axis point of the object (y4(z) = 0) in a direc- 
tipn determined by y3(z) = 1 (see fig. 1). 


6 (2g)= 


Yv(2Z,)' 


‘Fig. 1. Particular solutions of the integral equation. 


The situation of the image plane z = z, now follows from the 


M4 at c ; 
Uy—JAL hPL) (C—%) +f W(2) Efe, (64) (E—L1) (6,20) — 


z i , (10) 
1 — fade #(0)(C—e) +f (0) fl, PUC,) (C0) (Cpa) — 


cdrresponding to the edge of the field or by € = 00) 
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1. BREMMER 
85. The magnification. The magnification 47 can be derived from 
the solution v,(z) according to the formula: 


V1 e) 


Mx — 2 rigs ots) = Lyte). (1) 


2, (29) Vi(z,) 
The minus sign indicates that 17 represents a positive quantity if 
the electrons leaving an off-axis abject point in a meridional plane 
do intersect the axis once (or occasionally an odd number of times) 
before arriving at the image point (sec fig. 1). _ 

We could derive an expression for ,(z;) and thus for Af simply by 
substituting zs 2; in (8). The resulting formula, however, would . 
still depend on 2, which occurs both in the integrand and in the 
integration limits of the various terms. A more useful formula 
without having 2, in the integrands can be deduced from y,(z) by 
taking into account the Wronski property relative to y, and 4, viz. 

¥'1(2) ¥a(2) — ¥2(z) yi(z) = constant. 

The value of the constant follows from the substitution z = 2p, 
remembering that y,(%) = 1, y2(%) = 0 and y5(2)) = 1. Hence 

¥i(2) ya(2) — yal2) vi(2) = I. (12) 
A further application of (12) at z == 2, yields, in view of (11) and 
y2(2,) = 0, 
l 4(2,) 
Diese se aaa (13) 
M V4 (25) 
A series expansion for y3(z,) is casily derived from (9). The resulting 
final formula for 1/47 becomes 


1 Vz) 
Mo P49) 


{— I + f 4 BG) G29) = 


oe y 


— faz w7(e) | as, PE) EC — 8) G4) + 


70 0 
oa ¢ q 


+ f agH2C) f as, ACE) CE) | dE RENEE G2) —--f-04 


Ta 49 *0 


The magnification can be computed from this scrics without knowing 
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he exact situation of the image plane provided the ficld is negligible 
there. In fact, in that case the upper limit of integration 2; can be 
replaced once again by co or by the z-coordinate of the edge of the 
field. 

_ The expansion (13) can be transformed. into an other one, the 
terms of which show a more symmetrical structure. The new expan- 
sion is obtained by first shifting the factor of the integrand which 
contains z, to the front by means of a number of reversals of orders 
of integration. By additional reversals it proves to be possible to 
obtain integrals the upper limit of which equals z,. The final expan- 
sion then reads 


det Vi(z,) 
M ~~ V%(z5) 


ay 


{1 + dt, Pt) (i — 4) — 


ai 


—_ i de, 2(E) (C1 — 20) | di, (C2) (bo — 81) + 


f, 


+ fae, HG) (a fea) (ty fae, PEN o—t)—-. (15) 
20 a i 
We also give the corresponding scries in practical units for magne- 
tic lenses. For these lenses we have, according to (2), 


0.0221: -.5 
P(e) = = Pe) (16) 
if H and V are expressed in gauss and volt respectively. The first 
few terms of (15) thus yield for the magnification by magnetic lenses 

0.0221 f.. ; | 
V dé, H*(2,) (€; — %) — 


= 


0.000488 ae aye ie 
ayer © | de, HE) (61 20) / dE, H*(C5) (Co—b1) + a 


20 an 


[ase 420 fat? 0) (6) [ats HEN GLa» (17) 
49 : ot a 
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§ 6. The dependence of the osculating constants on the magnification. 
‘he function #?(z) of (2) decreases gradually on both sides of its 
maximum without being zero at the object and image planes, There- 
fore any clectron lens constitutes, strictly speaking, an immersion 
ens. As a consequence it is impossible to define focal and cardinal 
oints which are independent of the special position of these planes, 
hough the dependence is negligible provided 2?(z) is very small at 
P= zandz =z, Glaser) introduced so-called osculating focal 
und cardinal points which take into account the influence of h(z5) 
hnd &7(z,), His definitions are chosen such that the tangents at the 
ntersections of the trajectories with z = 4 and z = z, show the 
fame properties, in connection with the osculating elements, as the 
ectilinear trajectorics in ordinary optics. A disadvantage of Gla- 
;crs definitions is that they do not apply to all electron lenses. In 
what follows we start from slightly different definitions which are 
¢pplicable throughout and which are identical with Glaser’s 
dcfinitions.in so far as the latter exist. Our definitions correspond, ¢.g., 
fo that of the osculating elements used in planetary mechanics; in 
fhe latter these elements determine the cllipse which approximates 
est the actual trajectory in the vicinity of the point under consider- 
tion. 
The four osculating clements to be introduced here are the focal 
yoints F, and F, and the cardinal points Hy and H;, in object and 


Fig. 2. Showing the significance of the osculating elements. 


ifpage space respectively; the focal distances /, = 21g — 2p. and 
/]= Zr; — Zn; will be defined as positive if the situation corresponds 
tg that of fig. 2. The osculating clements are to be chosen in accord- 
apce with the concepts of elementary geomctrical lens theory, that 
is such that the image Q’ of the object point Q may be found as if 
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the electron trajectories were straight lines. The corresponding situa- 
tion as shown in fig. 2 recalls the well-known relations 
2H) — *ny = M (zr, — %o); 2; — &n, = Mer; — 2u1)- (18) 
In these equations 2, 2; and M are the only known quantities in our 
case. For a determination of the four unknown coordinates 2,,, 2m; 
ZH, -and zy, fixing the osculating elements we need two further 
relations. The latter follow from the condition that the same oscula-. 
ting elements should also apply to a second position of the object 
plane (and the corresponding image plane) that is infinitely close to 
the original one. This implies that (18) can be differentiated with 
respect to % while keeping 271,, 21; 2r, and zp, constant. In this way 
we get the additional relations 

dM dz; 
O= (zr, %0) dy ==; Pe = (24; — 2n;) 
By solving the four unknowns from the four equations (18) and (19) 
we obtain 


=% 1+ a 


dM 
dz 


(19) 


M dz; 
Sp, = 2, ; 
dM /dz, dzy 


M 
dM /dz,’ 
M+ M? " J4+M dz; 


(20) 
2n, = % + ar 


= 2;,-—- 


dM dz,’ on dM dz, day’ 


Further, the differéntial coefficient dz,/dz) can be connected with 
M as follows. We consider an object plane z= 2 + A% and a 
corresponding image plane z = z, + Az,, differing from the original’ 
ones (for which Azg= Az; = 0). The quotient of any two particular 
solutions of (1) assumes identical values at 2 = % + AZ and 
z= 42,-+ Az,, a property which expresses the homogeneous character 
of the magnication. An application of this property to the two 
special solutions y, and y, (see fig. 1) yields 


V1(% + Azo) volz; + 42%) — Yal%o + A%) Vi (2; + Az) = 0. 
Next we obtain from a differentiation with respect to 429 


Az, = __ MZ + Azo) ¥o(z% + A2;) — Val%o + Az) Vi (% + Az,) 
dA% Vi (Z le Azo) y9(2% + Az;) —Yol% + 4%) Vile + Az,)” 


We can pass.to the limit for Az, > 0, bearing in mind that y,(%) = 
== Yo(%) = 1 and yo(%) = Yo(z) = 0. Hence dz,/dzo = ¥1(2,)/¥o(2i), 
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or, in view of (I1) and (13), 
dz, F(z) 74 
: arf ie (21) 
(ls } (Zp) 
A substitution into (20) Icads to the following formulac if we intro- 
duce the focal distances /, = 2H, -—*r, and /, = zp; — %,, instead 
of zy, and Hy: 
ies MW F(z) \ AGB 
har 2p, = 2;-—- = - , 
aE ANE dc VI (eq) dM /dag 
ee MW? / F(z) \ M? 
o  dAL day’ "AN (e)5 dat fdz, 

These final expressions prove to be identical with those derived by 
Glaser. For a comparison with the formulae of this author, 
however, it is necessary to link JZ .-: —-9,(z,) to those particular solu- 
tions of (1) that correspond to electrons arriving at either side from 
infinity in the direction of the z-axis. 


$7. Explicit formulae for the osculating constants of magnetic lenses. 

I:xpansions similar to (15) can now be derived from (22) for the 
osculating elements. The evaluation can be based on (15) and its 
derivative with respect to zo. In the case of pure magnetic lenses we 
have F(z) = V'(z), so that the focal distances fy = /; = M?/(d.M/dzp) 
can be derived from 

l (147) a 

—_ dp 
Irom a differentiation of (17) we obtain the following expansion in 
practical apes 
l _ 0. 0221 


/ V 


fa: (2) — dz, (5) | dt, HP(2)(G—01) + 


/_ 1079.10 § 
ro. V3 


0.000488 
agi 


if as, 1) fazIPae—) fesse (53) 
70 q ti 


Poe (42) [es H?(&) (§ —2) (2; —¢) + 


bai) 


—IP(z,) M? 
ase Mea) fs FUN See eta) — (21) 

39 a 
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The second group of terms can be neglected if the field at the image 
plane is very small. Otherwise we may consider these terms as cor-- 
rection terms which can be taken into account by substituting appro- 
ximate values for M and z;. We emphasize that the number of inte-. 
grations occurring in the terms of the first group is lower than the 
corresponding number in the conventional formulae °). 

The positions of the focal and cardinal points can now be derived: 
from (17) and (21) with the aid of the following relations, which 
result from (18): 


Ge Sosy Sues eee ee 
Woo e Oe” ee 6 


The corresponding formulae for the elements of the image space 
are always found by interchanging the role of z and z;, which implies 
that M has to be replaced by 1/M. The explicit form of these for- 
mulae is much simplified if the additional terms arising from the 
field at the image (or object plane) may be neglected. In that case 
we obtain, e.g., the following expression for the position of the 
cardinal point in the object space: 


Zp, = % 1 


N 
4H, Fe 


in which | 


cA a4 


‘ 0.0221 : 
N= {at PQ) S——G— fac, HG) : At H(E5) (C2— 1) + 


29 29 fo 


0.000488 é ; /' 
° ye dg, HP(é,) c fac, ue eye) fate) eta ed 


a te 
| at, H(t) | at, H(by) (— ty) + 


% fi 
hi at ‘a 


0,000488 ee eee poe 
+f AE GY) AL, H (Eq) (Cg— 61) | Aly H7(E3)(Cn— 2) «+ = 
29 f bo 
For very fast electrons or very weak fields the first term of the 
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numerator and that of the denominator yield a sufficiently accurate 
approximation, viz. 


satire 
2uy™ Be, 


fas H(:) 


In this approximation both cardinal points have the same z-coor- 
dinate as the centre of gravity of the part of the curve representing 
H? as a function of z that is bounded by the curve, the z-axis and 
the ordinates z = z)andz =: 2, It is striking that in such expressions 
H? and not # constitutes the dominating quantity. 
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